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' Abstract 

o : 

t^J- , Non-uniform structures of the nucleon matter are expected at subnuclear densities 

and above the nuclear density: they are called nuclear pastas and kaon pastas, respec- 
tively. We numerically study these phases by means of the density functional theory 
\Ty . with relativistic mean-fields and the electric field; the electric field is properly taken 

into account. Our results demonstrate a particular role of the charge screening effects 
on these non-uniform structures. 
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1 Introduction 

There are expected various form of matter inside neutron stars (Fig. [TJ, many of which are 
associated with first order phase transitions (FOPT). Recently there are many studies of the 
mixed phases at these FOPT such as hadron-quark deconfinement transition PU 121 El E] , kaon 
condensation [5J El 13 El EE Ell HU U2] , color superconductivity [121 EH UH] , superfluidity in 
atomic traps [E], nuclear pastas [TTl [TBI [THl 1201 IZIl [221 1231 1231 , etc. Before the remark 
by Glendenning many authors used the Maxwell construction (MC) to get the equation of 
state (EOS) in thermodynamic equilibrium for FOPT. Nowadays there exists a view that 
not all Gibbs conditions can be satisfied in the description of the Maxwell construction in 
multi- component systems and the appearance of the mixed phases is inevitable, cf. P3 Ell- 
in this lecture we discuss two types of the mixed phases expected in nuclear matter; one is 
so called the "nuclear pastas" at subnuclear densities, and the other is the one following kaon 
condensation above the nuclear density. We shall see non-uniform geometrical structures in 
both cases. Consider a mixed-phase consisting of two phases in equilibrium, say I and II. The 
system should be totally neutral but composed of many particle species in chemical equilib- 
rium through weak processes. When we impose the baryon number and charge conservation 
on this system, we can easily see that only two chemical potentials are independent; one is 
baryon-number chemical potential \xb and the other is charge chemical potential \xq. Once 
\ib and hq are given, all the chemical potentials of particle species are determined. Then 
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each particle density may spatially change in the mixed phase by the strong and electroweak 
interactions, but chemical potentials must be constant over the whole space. The Gibbs 
conditions (GC) for thermodynamic equilibrium between two phases I and II are given by 



t ii pi = p u 



in this case, where T l and P % denote temperature and pressure of each phase, respectively. 

It has been claimed that 
usual MC satisfies only the first 
three conditions, and the fi- 
nal one is violated because it 
assumes the local charge neu- 
trality in each phase instead of 
the global charge neutrality. A 
naive application of the Gibbs 
conditions to separate bulk 
phases without the surface 
and the Coulomb interaction, 
demonstrates a broad density- 
region of the structured mixed 
phase [H Ej. When one takes 
into account the geometrical 
structures like droplet, rod and 
slab by extending the bulk cal- 
culation, one can see that the 
surface tension and Coulomb 
interaction determine their size 
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Figure 1: 
structure. 

[T71I21E1- However, the charge screening effect (caused by the rearrangement of the charged- 
particle distributions) should be very important when the typical size is of the order of the 
minimal Debye screening length in the problem. It may largely affect the stability condition 
of the geometrical structures in the mixed phases. We have been recently exploring the effect 
of the charge screening in the context of the various structured mixed phases [*H "it 112*1 l^7j. 

Our aim here is to investigate the non-uniform structures in nuclear matter numerically 
by means of the density functional theory with a relativistic mean field (RMF) model. Our 
framework allows one to determine the density profiles exactly without any sharp boundary 
and the surface tension used in the bulk calculations. It includes the Coulomb interaction 
in a proper way and we can fully take into account the charge screening effects. We shall 
figure out how the charge screening effects induce the rearrangement of the charged-particle 
distributions and thereby modify the results obtained by the bulk calculations. 



2 Nuclear pastas 

First we consider the nuclear pastas. At subnuclear densities, where pressure of the uniform 
nuclear matter is negative, so called "nuclear pastas" may appear [T7]. In view of the 
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phase transition, these can be regarded as the mixed phases following the liquid-gas phase 
transition in nuclear matter. 1 Stable nuclear shape may change from sphere to rod, slab, 
tube and to bubble with increase of the matter density. Pastas are eventually dissolved into 
the uniform matter (liquid phase) at a certain nucleon density below the saturation density, 
Po ~ 0.16 fm~ 3 . 

Existence of such pasta phases instead of the crystalline lattice of nuclei or the liquid npe 
phase would modify some important processes by changing the hydrodynamic properties and 
the neutrino opacity in the supernova matter and in the protoneutron stars. Also the pasta 
phases may influence neutron star quakes and pulsar glitches via the change of mechanical 
properties of the crust matter. 

A number of authors have discussed the nuclear pastas using various models ^7J QHl EH 
I2U1I2T1 1221 EMI l2o"| . but most of them have relied on the bulk calculations; roughly speaking, 
the favorable nuclear shape is determined by a balance between the surface and the Coulomb 
energies. Thus the Coulomb interaction as well as the surface tension has an crucial role in 
the non-uniform pasta structures. However, the treatment of the Coulomb interaction so far 
has been rather simple and the rearrangement effect on the density profiles of the charged 
particles by the Coulomb interaction has been discarded in the bulk calculations. The paper 
[26j discusses the effect of electron screening to demonstrate that it is of a minor importance, 
but the rearrangement of the proton density as a consequence of the Coulomb repulsion was 
not shown up in their study. 



2.1 Density functional theory with relativistic mean field 

Following the idea of the density functional theory (DFT) with the RMF model [SB] , we can 
formulate equations of motion to study non-uniform nuclear matter numerically. The RMF 
model with fields of mesons and baryons introduced in a Lorentz-invariant way is simple for 
numerical calculations, but realistic enough to reproduce the bulk properties of finite nuclei 
as well as saturation properties in nuclear matter. In our framework, the Coulomb interaction 
is properly included in equations of motion for nucleons, electrons, and meson mean-fields, 
and we solve the Poisson equation for the Coulomb potential Vcoui self-consistently with 
them. Thus the baryon and electron density profiles, as well as the meson mean-fields, are 
determined in a way fully consistent with the Coulomb potential. 

To begin with, we present the thermodynamic potential for the neutron, proton and 
electron system with chemical potentials /i n , \i v and /i e , respectively; 



n = n 



N 



M 



a 



(2) 



where 



n 



N 



a=p,n 



47T 3 



m* N + k 2 - p a v a 



with the local Fermi momenta, &F,a( r )( a — n iP)i f° r nucleons, 



(3) 



n 



M 



d 3 r 



(Va) 2 + m 2 a a 2 



+ U(a) 



1 Note that there is no uniform gas phase in nuclear matter at zero temperature. 



(4) 
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Table 1: Parameter set used in RMF in our calculation. 
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for electrons and the Coulomb potential, Vc ou i( r ), where v p (j) — pi p + Vc ou i( r ) — 9uiN^o( r ) — 
g pN R (r), z/ n (r) = p, n - g uN w (r) + fiVv^oO), m* N (r) = m N - g aN a(r), and the nonlinear 
potential for the scalar field, U(a) = ^brrij^^g^Na) 3 + |c(<7 CT ArO") 4 . Temperature T is kept to 
be zero in the present study. 

Here we used the local-density approximation (LDA) for nucleons and electrons. Strictly 
speaking, LDA is meaningful only if the typical length of the nucleon density variation is 
larger than the inter-nucleon distance. To go beyond LDA we must take into account some 
derivative terms with respect to particle densities, which can be easily incorporated in the 
quasi-classical manner by the derivative expansion within the density functional theory [2H] ■ 
In the case when we suppress derivative terms of nucleon densities they follow changes of 
the meson mean fields and the Coulomb field that have derivative terms. We must also 
bear in mind that for small structure sizes, quantum effects become prominent which we 
disregarded. Here we consider large-size pasta structures and simply discard the derivative 
terms, first-step calculation. 

Parameters of the RMF model are chosen to reproduce saturation properties of nuclear 
matter: the minimum energy per baryon —16.3 MeV at p = po = 0.153 fm~ 3 , the incom- 
pressibility K(p ) = 240 MeV, the effective nucleon mass m* N (po) = 0.78rriN', n^N = 938 
MeV, and the isospin-asymmetry coefficient a sym = 32.5 MeV. Coupling constants and 
meson masses used in our calculation are listed in Table 1. 

From the variational principle = (0* = a, Ro, u , V Co ui) and = (a = n,p, e), 
we get the coupled equations of motion for the mean-fields and the Coulomb potential, 



VV(r) = m>(r) + ^-^(pi s) (r)+pt s) (r)), (6) 

V 2 u; (r) = m^ofr) - g wN (p p {v) + p n (r)), (7) 

V 2 i? (r) = m 2 p Ro{r)-g pN {P P ( r )-Pn(r)), (8) 

Wcoui(r) = 47re 2 pch(r), (9) 

with the scalar densities p^(r),a = n,p, and the charge density, p c h(r) = Pp( r ) + Pe{ r )- 
Equations of motion for fermions yield the standard relations between the densities and 
chemical potentials, 



Pn = \Jkl (r)+m* N (r) + g wN ^o(r) - g pN Ro(r), (10) 
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P v = Pn - fie = \J ^F, P ( r ) + m v( r ) 2 + 5Wv^o(r) + g pN Ro{r) - Vcoui(r), (11) 
p e (r) = -(/i e - V C oui(r)) 3 /37T 2 , (12) 



where we have assumed the chemical equilibrium among nucleons and electrons. The baryon- 
number chemical potential /i# equals to p n and the charge chemical potential pq to p e . Note 
that first, the Poisson equation for the Coulomb field is a highly nonlinear equation in 
Vcoui( r ), since p c h( r ) in r.h.s. includes it in a complicated way. Secondly, the Coulomb 
potential always enters equations through the gauge invariant combinations p e — Vc ou i( r ) 
and p v + Vcoui( r )- Thirdly, solutions of these equations of motion attain the kinematical 
equilibrium, i.e. equal pressure in each spatial point, which is one of the Gibbs conditions. 



2.2 Numerical procedure 

To solve the above coupled equations numerically, we use the Wigner-Seitz cell approxima- 
tion: the whole space is divided into equivalent cells with a geometry. The geometrical shape 
of the cell changes: sphere in three dimensional (3D) calculation, cylinder in 2D and slab in 
ID, respectively. Each cell is globally charge-neutral and all the physical quantities in a cell 
are smoothly connected to those of the next cell with zero gradients at the boundary. Every 
point inside the cell is represented by the grid points (number of grids A^ gri d ^100) and the 
differential equations for fields are solved by the relaxation method under the constraints of 
the given baryon number and the global charge neutrality. 

To illustrate how to solve equations of motion for mean-fields, let us consider, for sim- 
plicity, two fields /i(r), f%{r) and their coupled Poisson-like equations under 3D calculation, 

V 2 /i(r) = mi 2 /i(r) + Wi[/i,/ 2 ], 

V 2 / 2 (r) = m 2 2 f 2 (r) + W 2 [f 1 ,f 2 ], (13) 

where Wi(i = 1,2) are functions of the fields fi and f 2 . Introducing a relaxation "time" t 
we solve 



dm 

dt 



(V^M-m^M-^)- (14) 



If coefficients q are appropriately chosen, the above /j will converge in time and we get the 
solution of Eq. (JT3J). 

Baryon densities are solved with the help of the "local chemical potentials" p a {f) (a — 
n,p), being different from the above introduced constant chemical potentials. Assuming 
/i n (r) being an increasing function of the neutron density p n (^) in Eq. (fTU|) . the relaxation 
equation for the neutron density, 



dt 



Cn(r) Pn(r)VVn(r), (15) 



is solved to equalize the local chemical potential fi n (r) in each point. The coefficient c n (r) 
is chosen to conserve the total neutron number. The proton density p p (r) is adjusted in the 
same way. When we impose the beta-equilibrium condition, proton and neutron densities are 
adjusted to achieve p n (r) = p p (r) + /i e (r). Finally we get densities p n (r) and p p (r) leading 
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Figure 2: Left: the density 
profiles of typical nuclei. The 
proton densities (solid curves) 
are compared with the exper- 
iment. Right: binding energy 
per nucleon and the proton ra- 
tio of finite nuclei. 



the constant chemical potentials p n (r) = and /x p (r) = \i v . Nevertheless, the basic idea is 
to obtain constant chemical potentials, [i a {f) = /i a (a = n,p) at the convergence. There is an 
exception: when there are some regions where p a {v) = 0, the local chemical potential /i a (r) 
is larger than the corresponding constant value in the regions with finite p a (r). 

The electron density p e {r) is calculated directly from Eq. ()12|). The value of /i e is adjusted 
at any time step to get global charge neutrality: we decrease \i e when total charge in a cell 
is positive and increase when it is negative. 

All the above relaxation procedures are performed simultaneously, which means that our 
numerical procedure always respect the Gibbs conditions for chemical potentials. 

3 Bulk properties of finite nuclei 

Before applying our framework to the problem of the nuclear pastas , we check how it works 
to describe the bulk properties of finite nuclei. In this calculation for simplicity we assume 
the spherical shape of nucleus. The electron density is set to be zero. Therefore the global 
charge neutrality condition is not imposed. 

In Fig. El (left panel) we show the density profiles of some typical nuclei. One can see how 
well our framework may reproduce density profiles. To get a better fit, especially near the 
surface, we could include the derivative terms of the nucleon densities, as we have mentioned. 
Fine structures seen in the empirical density profiles, which come from the shell effects (see, 
e.g., a proton density dip at the center of a light 16 nucleus), cannot be reproduced by the 
mean-field approach. The effect of the rearrangement of the proton density can be seen in 
heavy nuclei; protons repel each other that enhances their concentration near the surface of 
the heavy nuclei. This effect is analogous to the charge screening effect in a sense that the 
proton distribution is now changed not on the scale of the radius of the nucleus, as for bare 
Coulomb field, but on another length scale, that we will call the proton Debye screening 
length, see Eq. (|T7|) below. It has important consequences for the pasta structures since 
typically the proton Debye length is less than the droplet size. The stable value of the 
proton fraction Y p = Z/A [Z and A are proton and total baryon numbers, respectively) is 
obtained by imposing the beta equilibrium condition fi n = p p for a given baryon number. 
Figure 121 (right panel) shows the baryon number dependencies of the binding energy per 
baryon and the proton ratio. We can see that the bulk properties of finite nuclei (density, 
binding energy, and proton to baryon ratio) are satisfactorily reproduced for our present 
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Figure 3: Examples of the density profiles in the cell for symmetric nuclear matter with 
1^=0.5 (left) and for asymmetric matter with Y p = 0.3 (center) and 0.1 (right). 



purpose. 

Note that we must use a slightly smaller value of the sigma mass i.e. 400 MeV, than that 
one usually uses to get an appropriate fit. If we used a popular value m a ~ 500 MeV, finite 
nuclei would be overbound by about 3 MeV /A. The actual value of the sigma mass (as 
well as the omega and rho masses) has little relevance for the case of infinite nuclear matter, 
since it enters the thermodynamic potential only in the combination C a = gaN/^a- However 
meson masses are important characteristics of finite nuclei and of other non-uniform nucleon 
systems, like those in pastas. The effective meson mass characterizes the typical scale for 
the spatial change of the meson field and consequently it affects, e.g., the surface property 
of the given nucleon structure and it influences the value of the effective surface tension. 



4 Non-uniform structures in nuclear matter 
4.1 Nucleon matter at fixed proton fraction 

First, we focus on the discussion of the behavior of the nuclear matter at a fixed proton 
fraction Y p . Particularly, we explore proton fractions Y p = 0.1, 0.3, and 0.5. The cases 
Y p = 0.3 - 0.5 might be relevant for the supernova matter and for newly born neutron stars. 
2 Figure El shows some density profiles inside the Wigner-Seitz cells as functions of the 
radial distance from the center of the cell. The geometrical dimension of the cell is denoted 

2 Note that we have no relation among nucleon chemical potentials and electron and three chemical 
potentials become independent in this case, different from the usual (3 stable matter. Accordingly we impose 
one more condition on the system as a Gibbs condition. 
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Figure 4: Binding energies per nucleon and the cell size for symmetric nuclear matter with 
Y p =0.5 (left), and for asymmetric matter with Y p = 0.3 (middle) and 0.1 (right). 



as "3D" (three dimensional), etc. The cell boundary is indicated by the hatch. From the 
top to the bottom the configuration changes like droplet (3D), rod (2D), slab (ID), tube 
(2D), and bubble (3D). Thus the nuclear pastas are clearly manifested. For the lowest Y p 
case (Y p = 0.1), the neutron density is finite at any point: the space is filled by dripped 
neutrons. For a higher Y p , the neutron density drops to zero outside the structure. The 
proton density always drops to zero outside the nuclear lump. The charge screening effects 
are pronounced, especially due to protons. Protons repel each other and thereby the proton 
density profile substantially deviates from the step-function: the proton number is enhanced 
near the surface of the nuclear lump. The electron density also becomes non-uniform by the 
rearrangement effect. This non- uniformity of the electron distribution is more pronounced 
for a higher Y p and a higher density. 

EOS of the nuclear pastas is shown as a function of the averaged density in Fig. 0] (upper 
panels). The energy E — itlb includes the kinetic energy of electrons, which makes the total 
pressure positive. The lowest-energy configuration is selected among various geometrical 
structures for a given averaged density. The most favorable configuration changes from the 
droplet to rod, slab, tube, bubble, and to the uniform one (the dotted thin curve) with 
increase of the density. We can see that the appearance of the nuclear pastas results in a 
softening of the EOS: the energy per baryon gets lower up to about 15 MeV/A compared 
to the uniform case. The lower panels in Fig. |U show the cell sizes i? ce ii and structure sizes 
i?d as functions of the averaged density. The size is defined here by way of a density 
fluctuation as 



Rd 



R ce w p , (for droplet, rod, and slab) 

(pi) 



2 x (16) 



R ce w ^1 — ~j^y j 5 (f° r tube and bubble) 



where the bracket "()" indicates the average inside the cell. Dashed curves show the Debye 
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screening lengths of electron and proton calculated as 



A 



(e) 
D 



-Aire 2 




-1/2 



A 



(p) 

D 



Aire 



dfip 



-1/2 



(17) 



respectively, where p^ v is the proton density averaged over the nuclear lump and p^ v is the 
electron density averaged inside the cell. Note that these quantities are obviously gauge 
invariant. Numerically, the cell size R ce n for droplet, rod, and slab configurations at Y p = 0.5 
and 0.3 are shown to be close to the Debye screening length of electron. For Y p = 0.1, in all 
cases _R ce ii is substantially smaller than A^ and thereby the electron screening should be much 
weaker. In all cases, except for bubbles (at Y p = 0.5 and 0.3), the structure size Rj are smaller 
than A^. This means that the Debye screening effect of electrons inside these structures 
should not be pronounced. For bubbles at Y p = 0.5 and 0.3, A^, is substantially smaller 
than the cell size and the electron screening should be significant. For Y p — 0.5, 0.3, 0.1 in all 
cases (with the only exception Y p = 0.1 for slabs), the value A^ is shorter than R^, which 
means that the rearrangement of proton density is essential for the structures of the nuclear 
pastas, as it is indeed seen from the Fig. 01 



Problem: Evaluate the electron Debye screening length in the case of p a e 
by using the above expression Eq. (fTYj) for massless electrons. 



-0.5p , 



Using the baryon density and the structure 
size from Fig. 01 one may estimate the atomic 
number of the structure. In the case of droplets 
and for Y p = 0.5 the atomic number of the 
droplet is ~ 25 at lower density limit and ~ 65 
at the maximum density of the droplet phase 
p^f ~ 0.025 far 3 . 

4.2 Nuclear matter in beta equi- 
librium 

Next, we explore the nuclear pastas in beta 
equilibrium. The droplet structure appears, 
which is quite similar to the above considered 
case of the fixed proton ratio Y p = 0.1. The 
apparently different feature in this case is that 
only the droplet configuration appears as a non- 
uniform structure. It should be noticed, how- 
ever, that the presence or absence of the each 
pasta structure may sensitively depend on the 
choice of the effective interaction. 

In Fig. El we plot EOS (top), the structure 
size (middle), and the proton ratio (bottom). 
The difference between EOS of uniform matter 



full calculation 




Figure 5: Binding energy (top), structure 
size (middle), and proton ratio (bottom) 
in the cell for nuclear matter in beta equi- 
librium. 
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Figure 6: Comparison of the density profiles for different treatments of the Coulomb inter- 
action. From the left: "full" calculation, "no Coulomb" calculation, and "bulk" calculation. 
The proton ratio for all cases is Y p =0.5. 



and that of non-uniform one is small, while the proton ratio is significantly affected by the 
presence of the pasta at lower densities. The droplet radius and the cell radius in the middle 
panel of Fig. 03 are always smaller than the electron Debye length A# , and thereby the effect 
of the electron charge screening is small. On the other hand, the proton Debye length A^ 
is comparable with the droplet radius at all densities, which demonstrates the relevance of 
the proton screening. 



5 Charge screening effect in nuclear pastas 

In this section we explore the effect of the charge screening on the nuclear pastas. Here we 
focus on the matter with fixed proton ratio Y p = 0.5 since the Coulomb effects should be most 
pronounced in this case. We compare three kinds of calculations with different treatments 
of the Coulomb interaction. One is the "full calculation" which we have presented above. 
The second calculation ("no Coulomb") is performed by totally discarding the Coulomb 
potential Vc ou i in equations of motion. After getting the density profiles, the Coulomb 
energy, being evaluated using charge densities thus determined, is added to the total energy. 
Note that this calculation is similar but not the same as the bulk calculations: in the latter 
a sharp boundary and the surface tension are introduced, and the particle density profile is 
assumed to be constant in each phase. The third is the bulk calculation, in which we use a 
completely uniform electron density, the baryon density distribution of a step function, and 
the surface tension introduced by hand. To determine the geometrical structure we have 
used the equilibrium condition of baryon chemical potential and the pressure calculated the 
present RMF model, and the relation between surface energy E snT f and the Coulomb energy 
-Ecoui [17J as 

Egurf = 2i?Coul- (18) 
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p/po u=p/p 

Figure 7: Onset mechanism of kaon conden- Figure 8: Original EOS and the Maxwell 
sation. construction. 

We have used the surface tension parameter r = 1.03 MeV/fm 2 to fit the liquid-drop binding 
energies of finite nuclei. 

EOS as a whole (upper panels in Fig. El) shows almost no dependence on the calculation. 
This agrees with a general statement that the variational functional is always less sensitive 
to the choice of the trial functions than the quantities depending on these trial functions 
linearly. 

Nevertheless, the density region of pasta structures and the sizes of the structures (lower 
panels in Fig. [BJ) especially for tube and bubbles are different. In fact we don't observe tube 
and bubble configurations in the present bulk calculation. In Ref. and others, however, 
they have reported the appearance of full pasta structures by the bulk calculation. The 
most important origin of this discrepancy may come from the difference in the nuclear EOS 
used in the calculation. Also the surface tension is crucial; if we use smaller value of r, e.g. 
0.3 MeV/fm 2 , we observe full pasta structures appear in wider density range. If we take 
r = the mixed phase spreads from zero to the saturation density po without any specific 
geometry. From the bulk calculation we see that the surface tension plays a crucial role in 
the appearance of pasta structures. 

Comparing the "no Coulomb" calculation with "full calculation", we don't see large 
difference. However, precisely looking, the density region of pasta in "full calculation" is 
slightly larger. The density dependence of the cell size is also different. In the case of "no 
Coulomb" all the pieces of lines for i? ce n are monotonously decreasing. The behavior of R ce \\ 
in the full calculation is rather similar to the bulk calculation. 



6 Kaon condensation in high-density matter 

Next we explore the high-density nuclear matter in beta-equilibrium, which is expected to 
exist in the inner core of neutron stars. Kaons are the lightest mesons with strangeness, and 
their effective energy is much reduced by the kaon-nucleon interaction in nuclear medium. 
For low-energy kaons the s-wave interaction is dominant and attractive in the 1=1 channel, 
so that negatively charged kaons appear in the neutron-rich matter once the process n — > 
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Table 2: Additional parameters used in RMF with kaon. The kaon optical Uk is defined 
by U K = goKO + fiW^o- 



M~A)[MeV] m K [MeV] 9u]K g pK U K (p ) [MeV] 
93 494 g pN -120 - -130 



p + K becomes energetically allowed. Since kaons are bosons, it causes the Bose-Einstein 
condensation at momentum p = [2H|. The threshold condition then reads 

e_(p = 0) = p n - p p = p e , (19) 

which means the kaon distribution function diverges at p = (Fig. [TJ). 

If kaon condensation occurs in nuclear matter, it has many implications on compact 
stars; softening of EOS may give the possibility of the delayed collapse of a supernova to the 
low-mass black hole, and the nucleon Urea process under background kaons may give a fast 
cooling mechanism of neutron stars |3UJ ED E2 

Since many studies have shown that kaon condensation is of the first order, we must 
carefully treat the phase change (Fig. |HJ). In the following we discuss the kaon mixed phase 
in a similar way to nuclear pastas. There have been some studies about the mixed phase in 
the kaon condensation El 13 El El EEE HH H2j • In ref . [TTj the charge screening effect has 
been also studied, but all the equations of motion have not been solved self-consistently. 

To incorporate kaons into our calculation, the thermodynamic potential of Eq. (J2J) is 
modified as 



n 



tt = tt N + tt M + tt e + tt K , (20) 
i< I u 1 — -ii'k ••• U'i< - » c.ui ••• /aw^'h ••• //„/v/»o.i -1- — — 7 (21) 



f d 3 r 






to 



where m* K = rax — gaK^, Pk = p e , and the kaon field K = /kO/V^ (fx ■ Kaon decay 
constant). 3 The equations of motion are then similar to Eqs. (5) - (11) given for nuclear 
pastas except kaon contributions • Additional parameters concerning kaons are presented 
in Table 2. 

If Glendenning's claim were correct, the structured mixed phase would develop in a 
broad density range from well below to well above the critical density |S1E1IZ|- In this density 
interval the matter should exhibit the structure change similar to the nuclear pastas [27J: the 
kaonic droplet, rod, slab, tube, bubble. Actually we observe such structures (kaonic pastas) 
in our calculation. In the top (a) and the middle (b) panels of Fig. Elwe show EOS; pieces of 
solid curves indicate the energetically most favored structures, while the dotted curve EOS 
of the uniform matter. One can see the softening of EOS by the appearance of kaonic pastas. 
In the bottom panel (c) plotted are the size of the kaonic lump or hole and the cell size. 



3 We here consider a linearized KN Lagrangian for simplicity, which is not chiral-symmetric. 
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The dashed lines show the Debye screening 
lengths of the electron and the proton, A^ and „ 

, respectively. In most cases A^ is less than 
the cell size R ce u but it is larger than the struc- 
ture size R. However the proton Debye length 
\p is always shorter than R ce \\ and R. When 
the minimal value of the Debye length inside 
the structure is shorter or of the order of R, the 
charge screening effects should be pronounced 
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Figure 9: Top: binding energy per nu- 
cleon of nuclear matter in beta equilib- 
rium. Middle: pressure. Bottom: struc- 
ture size R (thick curves below) and the 
cell size _R ce ii (thick curves above). 



7 Charge screening effect 
in the kaon mixed phase 

To demonstrate the charge screening effect on 
the kaonic mixed phase and discuss differences 
among various treatments , we show pressure 
P (= -n/V) in the middle panel (b) of Fig. M 
and the phase diagrams in the fiB-^e plane in 
Fig. EH From Fig. El (b) we can clearly see 
that our results give the similar pressure to the 
one given by MC, while the bulk calculation, 
where no Coulomb interaction or surface energy 
is taken into account, gives a wide density re- 
gion for the mixed phase. 

We can get more insight about the role of 
the Coulomb interaction in Fig. fTOJ left panel in Fig. [TO] exhibits the full calculation, while 
in the right panel we show the case, when the electric potential is discarded in determining 
the density profile ("No Coulomb") and the Coulomb energy, using the density profile thus 
determined, is then added to the total energy. We use "No Coulomb" in the same meaning 
as in the nuclear pastas. We see that in "No Coulomb" case the pieces of solid curves lie 
between two curves given by the bulk calculation (indicated by "Gibbs"), where the Gibbs 
conditions are imposed disregarding the surface tension and the Coulomb interaction, and 
by the Maxwell construction (indicated by "Maxwell"). The "Full calculation" case is more 
close to the one given by the Maxwell construction. As follows from the density profiles 
shown in Fig. the local charge neutrality is more pronounced in the case of the "Full 
calculation" (smaller difference of kaon and proton densities). These results suggest that the 
Maxwell construction is effectively meaningful owing to the charge screening effects. 



8 Summary and concluding remarks 

We have discussed two kinds of the non-uniform structures in nuclear matter, nuclear pastas 
at subnuclear densities and the kaon mixed phase above the nuclear density, which may arise 
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M-b [MeV] ji B [MeV] 

Figure 10: Phase diagram in the /ie-/i e plane, where \ib and /i e denote baryon- and charge- 
chemical potentials, respectively. Left: full calculation. Right: the electric potential is 
discarded in determining the matter structure. Curves calculated using Gibbs conditions 
disregarding finite size effects and the Coulomb interaction effect ( "Gibbs" ) and that for the 
Maxwell construction ("Maxwell") are also presented for comparison. 



Full calculation ' 'No Coulomb' ' 




Figure 11: Comparison 
of density profiles of kaonic 
matter calculated with and 
without the Coulomb inter- 
action. The cell size is not 
optimized in this compari- 
son. 



r [fm] r [fm] 



as consequences of FOPT with many particle species, and elucidated the charge screening 
effect. Using a self-consistent framework based on density functional theory and relativistic 
mean fields, we took into account the Coulomb interaction in a proper way and numerically 
solved coupled equations of motion to extract the density profiles of nucleons, electrons and 
kaons. 

First of all, we have checked how realistic our framework is by studying the bulk properties 
of finite nuclei, as well as the saturation properties of nuclear matter, and found that it can 
describe both features satisfactorily. One could still improve the consideration by including 
the gradient terms of the nucleon densities, which may give a better description of the nuclear 
surface. 

In isospin-asymmetric nuclear matter with fixed proton to atomic number ratios, we have 
observed the "nuclear pastas" with various geometrical structures at sub-nuclear densities. 
These cases are relevant for the discussion of the supernova explosions and for the description 
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of the newly born neutron stars. The appearance of the pasta structures significantly lowers 
the energy, i.e. softens the equation of state, while the energy differences between various 
geometrical structures are rather small. 

By comparing different treatments of the Coulomb interaction, we have seen that the 
self-consistent inclusion of the Coulomb interaction changes the phase diagram. In partic- 
ular the region of pasta structure is broader for "full" calculation compared to that with 
simplified treatments of the Coulomb interaction which have been used in the previous stud- 
ies. The effect of the rearrangement of the proton distributions on the structures is much 
more pronounced compared to the effect of the electron charge screening. The influence of 
the charge screening on the equation of state, on the other hand, was found to be small. 

We have also studied the structure of the nucleon matter in the beta equilibrium. We 
have found that only one type of structures is realized: proton-enriched droplets embedded 
in the neutron sea. No other geometrical structures like rod, slab, etc. appeared. 

We have discussed how the geometrical structures manifest in the context of kaon con- 
densation. Our framework can be extended to include kaons straightforwardly. We have 
discussed the effect of the charge screening in this case. Since the kaon mixed-phase appear 
at high-densities, we see that changes are more remarkable than for the "nuclear pastas" at 
subnuclear densities |27] . The density range of the structured mixed phase is largely limited 
by the charge screening and thereby the phase diagram becomes similar to that given by the 
Maxwell construction. Although the importance of such a treatment has been demonstrated 
for the hadron-quark matter transition [3J H], one of our new findings here is that we can 
figure out the role of the charge screening effect without introducing an "artificial" input of 
the surface tension. In our study we have used the one-boson-exchange interaction for the 
KN interaction. On the other hand the bulk calculation, where no Coulomb interaction nor 
surface tension is included, can not give the mixed phase in the chiral model [HI E] • It should 
be interesting to study the effects of the charge screening to see whether the chiral model is 
not thermodynamically well-defined. 

In application to the newly formed neutron stars like in supernova explosions, finite 
temperature and neutrino trapping effects become important, as well as the dynamics of 
the first order phase transition with formation of the structures. It would be interesting to 
extend our framework to include these effects. 
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